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Abstract: Triple gauge boson associated production at the LHC serves as an interesting 
Q-i channel to test the robustness of the Standard Model. Any deviation from its SM prediction 

may indicate possible existence of relevant new physics, e.g., anomalous quartic gauge 
boson couplings. In this paper, a Monte-Carlo feasibility study of measuring W + W~^f 
production with pure leptonic decays and probing anomalous quartic gauge-boson (e.g., 
) couplings, is presented in detail for the first time, with parton shower and detector 
simulation effects taken into account. Our results show that at the y/s = 14 TeV LHC with 
an integrated luminosity of 100 (30) fb _1 , one can reach a significance of 9 (5) a to observe 
the SM W^ + W^~7 production, and can constrain at the 95% CL the anomalous WW'yy 
coupling parameters, e.g., Oq^/A 2 (see Ref. [14] for their definitions), below 1 x 10 _5 GeV" 2 , 
respectively. 
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1 Introduction 

The Standard Model (SM) of particle physics has reached a great success below the TeV 
energy scale, especially after the discovery of the 125-126 GeV Higgs-like boson this year [1- 
4]. So far, we haven't yet observed any significant deviation from the SM. However, there 
are many reasons to expect new physics beyond the SM appearing at the LHC (TeV) energy 
scale, such as the demand for dark matter candidates and the quest to understand large 
hierarchy between the electroweak and Planck scales. Further confirmation of the SM or 
uncovering new physics beyond the SM thus become as urgent goals for both theorists and 
experimentalists at the LHC era. 

Within the framework of the SM, gauge boson self-interactions are completely deter- 
mined by the SU(2)l ® U(l)y gauge symmetry, thus direct investigation of gauge bo- 
son self-interactions provides a crucial test on the gauge structure of the SM. Moreover, 
since the longitudinal components of and Z° result from the spontaneously symmetry 
breaking, these kinds of study may also be important to explore the electroweak symmetry 
breaking (EWSB) mechanisms. 

For the triple gauge boson associated production which is of our interest in this paper, 
extra contributions other than the SM predictions can be induced by possible new physics, 
which can be expressed in a model independent way by high-dimensional operators which 
lead to anomalous triple or quartic gauge boson couplings (TGCs or QGCs). Compared 
with the gauge boson pair production channel for TGCs measurement, triple gauge boson 
production, although suffered from lower cross sections and more complicated final state 
topology, is crucial for testing QGCs, especially in the cases as discussed in e.g., Ref. [5, 6], 
that it is possible the QGCs deviate from the SM prediction while the TGCs do not, 
assuming e.g., the exchange of extra heavy boson generates tree-level contributions to four 
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gauge boson couplings while the effect on the triple gauge vertex appears only at one-loop 
and is consequently suppressed. 

While considerable efforts have been made on probing TGCs, less are for QGCs. Pre- 
vious Monte-Carlo (MC) and experimental studies on QGCs have been done at ej and 77 
colliders [7, 8], linear colliders [5, 9-14], and hadron colliders [6, 15-19]. Direct constraints 
on QGCs are available from the LEP collider searches via W + W~"f [11, 13], 77^(ij) [14] 
and Z77 [12] channels, e.g., constraints at the 95% CL on WW~f~f QGC parameters a^/A 2 
are at the order of 10 _2 GeV" 2 [14], where (Zq^, £1X6 dimensionless coupling constants and the 
A stands for a new physics scale. 

With the excellent performance of the LHC and its upgrade plan to higher collision 
energy or luminosity in the next few years, we expect that the LHC may play a crucial role 
in determining the QGCs, as shown in Refs. [6, 15, 16] through MC simulation studies on 
direct VF77 , VBF 77, Z7 and WW channels, from which constraints on a^/A 2 can be 
lowered to reach the order of 10 -5 — 10 _6 GeV~ 2 with integrated luminosity of 30 — 100 fb _1 . 

In this paper we focus on the MC feasibility study at the LHC of measuring W + W~^f 
production with pure leptonic decays, and probing the anomalous quartic gauge-boson 
coupling (WW'yy as an example), with parton shower and detector simulation effects 
taken into account. We begin by specifying the effective Lagrangian related to the QGCs 
in Sec. 2, and then describe the framework of our simulation studies in Sec. 3, with the 
selection cuts demonstrated in Sec. 4. We present the LHC sensitivities on W + W~^f 
production with W + W~ leptonic decays and the QGCs subsequently in Sec. 5. Finally we 
conclude in Sec. 6. 

2 Effective Interactions for Photonic QGCs 

An effective Lagrangian can be constructed in a model independent way for the anomalous 
quartic couplings, assuming that new Physics beyond the SM keeps SU(2)l®U(1)y gauge 
invariance. The Lagrangian can usually be written down either linearly or non-linearly [5, 
10]. The lowest order genuine QGC operators are dim-6 for the non-linear representation 
and dim-8 for the linear representation [6, 10]. 

Although it is now more preferable to work in the linear context due to the recent 
discovery of Higgs-like boson, we note here that the two ways indeed involve similar QGC 
operators [6, 10]. For example, there are only two basic Lorentz structures for the WW'yy 
QGCs which are of our interest in this paper: 



Wo = ~F^F^W +a W~, (2.1) 



W] = -^F^F^{W +V W- + W~ V W+), (2.2) 
Accordingly, the effective interactions can be expressed by the above operators as [10] 
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C = k^W^ + k^W2. (2.3) 

Here k^ c are dimensional QGC parameters which can be compared with Oq^/A 2 in e.g., 
Refs. [14, 16]. 

One should also note that the effective Lagrangian leads to tree-level unitarity viola- 
tion. To avoid that, we take the following common used dipole form factor 

k w 

(l + 5 °/A2)2' (2 ' 4) 

where s is the triple gauge boson invariant mass and the parameter A is fixed to be 2 TeV. 

Finally, we want to point out that the linear way in principle leads to correlations 
among different QGC operators and thus can be more complicated than the non-linear 
way. However, for simplicity and convenience of comparison with previous experimental 
results, we still assume only photonic WW77 QGCs and omit all others, from a bottom-up 
point of view. 



3 Event Simulation 

The characteristic signal we are interested in contains two well identified lepton (electron 
e, or muon fi) in association with large missing transverse energy p T . In Fig. 1, we show 
examples of Feynman diagrams for the W + W~j productions at the LHC in the di-leptontic 
final state lulwy, with I = e, \i and r. Note r decays into e, /i at the ratio of about 35% and 
is handled with TAUOLA [21]. Five main background processes are considered: Z7, ZZ7, 
ZW'y, tW"f and tt'y, where ti'y is the dominant one. Here we don't consider backgrounds 
with photons from jet fragmentation, in which the photons tend to be close to jets and the 
contributions can be suppressed efficiently via photon isolation cuts (see e.g. Ref. [20]). 

Attention should be paid to Fig. 1(c), as it can also be seen as the initial and final 
state radiations (ISR and FSR) from the WW production process, generated by Pythia. 
However, the ISR and FSR approximations in Pythia should break down for hard or wide 
scattering photon, e.g., when the transverse momentum of 7, Pt~, is large. Note also this 
subset of contributions to VF + VF~7 is not related to QGCs, thus it would be interesting 
and important to show the overall W + W~~f results subtracting the contributions of the 
ISR/FSR approximations of Fig. 1(c), which we denote as pure_Vs: 

pure_Vs = WW7-ISR/FSRWW. (3.1) 

Our work is mainly carried out under MadGraph/MadEvent [22, 23, 27]. The 
photonic QGC effective Lagrangian 2.3 has been implemented into MadGraph using 
the FeynRules [24]-UFO [25]-ALOHA [26] framework. Signal and background events are 
then generated with MadGraph [22, 23] and MadEvent [27], processed subsequently by 
Pythia 6 [28] for parton showering and hadronization. Finally, the events are passed to 
Delphes version 2.0 [29] for detector simulation, where we focus on the CMS detector at 
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Figure 1. Example Fcynman diagrams contributing to W + W 7 productions at the LHC 

the LHC. Finally, the analysis is performed with the program package ExRoot Analysis [30] 
and ROOT [31]. The work flow has also been used in our previous study on semi-leptonic 
decayed WW simulations [32] . 

4 Event Selection 

We choose the following pre-selection cuts to generate unweighted events at parton level 
with MadGraph/MadEvent to interface later with Pythia and Delphes , 



for the signals and backgrounds listed in Sec. 3, where 77 is the pseudo-rapidity and 4> is the 
azimuthal angle around the beam direction. Note, however, for the backgrounds ZZ'y and 
ZW^f where one lepton is misidentified, we don't require any of the above cuts on leptons 
in order not to make bias. 

Moreover, in the hard process generation with MadGraph/MadEvent we adopt 
the CTEQ6L1 parton distribution functions (PDFs) [33] and set the renormalization and 
factorization scales as the transverse mass of the core process. 

Tighter cuts are then imposed on the reconstructed objects in the Delphes settings 
cards, 

• Pr e ,M > 20GeV, and |r? e;/vr | < 2.4. 

• Jets are clustered according to the anti— kt algorithm with a cone radius AR = 0.5. 
Moreover, Ptj > Pj^j (25 GeV by default) and \r]j\ < 5 are required. 

Other high level cuts are set in the analysis step as following: 

• (1) One and only one photon is allowed per event, and P Tl > Pf * with = 





• P Tl > 20 GeV, and |r/ 7 | < 2.5, 



• R^i = ^Ar^ + A^ > 0.4, 



40 GeV by default. 
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• (2) Two and only two leptons with opposite sign of charge. 

• (3) $ T > 40 GeV. 

• (4) Rji and R~j are larger than 0.5. 

• (5) To suppress top quark related backgrounds, we exclude events with b-tagged jet 
(with default Delphes's setting, e.g. b-tagging efficiency as 40%). 

• (6) To suppress backgrounds with Z boson leptonic decay, we require \mu — Mz\ > 
10 GeV. 

In Delphes, photons and charged leptons may overlap with the jet collections: Delphes 
first reconstructs photons and leptons based on MC information, and then jets which can 
be seeded from the already reconstructed photon or leptons. In our analysis, we clean the 
lepton collections from jets by requiring the Delphes's calculated "EhadOverEem" (the 
energy deposition in the Hadron Calorimeter over the one in the Electromagnetic Calorime- 
ter) smaller than 1. Moreover, we removed any jet which has Rj~ < 0.001 as it would be 
indeed most like a photon. 

5 Numerical Results 
5.1 W + W~y production 

As a first step, we are interested in estimating the feasibility of observing triple gauge boson 
W + W~"f productions at the LHC, before going into anomalous QGCs. As mentioned 
before, we are also interested in comparing overall W+W'-y results with the ISR/FSR 
ones from WW processes (see the VFVl^-column in Table 1). 

To optimize our results, we introduce further the following 3 requirements individually, 
in addition to all the cuts mentioned in Sec. 4: 

• (A) Veto events with jets of which P?j > P?j> thus we have P™j < Pt j < P^j: 

• (B) Vary P™%, 

• (C) Vary Pj."- and require jet numbers rij to be or 1. 

We are setting the special cuts of (A) and (C), as we would like to suppress more 
top-quark related background while keeping high signal efficiency. Compared with signal, 
ttj and tW"f are different in two sides: they are QCD processes and tend to radiate more 
jets; the hard physics scale is higher and thus one or more jets can be harder than the jets 
in signal events. 

We list the event numbers for the signal and backgrounds in Table 1, with the op- 
timized parameters (optimized for pure_Vs contributions) for the above 3 cases: (A*) 
pUP. = 60 GeV, (B*) P Tl > 80 GeV, and (C*) P Tj > 25 GeV with nj < 1. Related 
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K-factors for the signal and backgrounds are also listed with references in Table 1. Cor- 
respondingly, the significances are shown in Fig. 2, calculated with Eq. (5.1) [34, 35]. 



Signif = ^2ln{Q), Q = (1 + N s /JV& ) Nobs exp(—N s ), (5.1) 

From Table 1), we can see that tij is the dominant background. It is also interesting 
to notice the QED ISR/FSR contributions from WW get decreased a lot when P™* is set 
to a high value as 80 GeV, as the QED radiation approximation in Pythia breaks down 
for hard photons. 



Processes 


Cross section 


K-factor 


Events 


[fb] 


[Ref.] 


(A* ) P%* = 60 GcV 


(B*) = 80 GcV 


(C*) rij = 0, 1, P Tj > 25 GcV 


W + W~~{ 


18.286 


2.0 [36] 


95.818 


58.880 


114.84 


I(F)SR WW 


3114.1 


1.5 [37] 


35.812 


4.6712 


54.498 


Z 1 


4107.2 


1.5 [38] 


61.608 


47.232 


57.501 


ZZ^ 


45.818 


1.3 [36] 


0.2779 


0.1985 


0.2780 


W ± Z-/ 


1.3698 


1.5 [39] 


0.8903 


0.5739 


1.0068 


tij 


170.22 


1.9 [40] 


88.830 


73.738 


60.801 


tw^ 


26.858 


1.0 [37] 


17.905 


11.442 


16.527 



Table 1. Cut flow at the LHC with ^/s = 14 TeV and integrated luminosity of 100 fb 



More details can also be checked in Fig. 2, showing the significance dependences on (A) 
jet vetoing cut P^p (B) P^, and (C) P^j. Note we also give the pure_Vs-curves to show 
the results after subtracting ISR/FSR contributions from WW processes, as mentioned 
above. 

(A) , ti'j and tWj backgrounds tend to have harder jet than the signal. Increasing the 
Jet vetoing cut P^, p - from low value as 30 GeV, at first can veto more top background while 
less signal events, and thus enhance the significance for both VF + VF~7 and pure_Vs cases. 
However, the gain gets small when P^ p - is as large as 60 to 100 GeV (as top decayed jets 
are now soft compared with Py?), for pure_Vs and W + W _ 7 , respectively. 

(B) . Increasing Pr^ lower cut decreases the overall significance of W + W _ 7 , as the 
events get decreased for both signal and backgrounds. However, the pure_Vs-significance 
increased, as the percentage of pure_Vs over W + W~j is enhanced. 

(C) . Increasing jet reconstructing cut Pj>"-, increases the 0-jet contributions while 
decreases the 1-jet ones, as expected. The overall 0+1 jet significances are decreased, as 
less ti'y and tWj events (which have more QCD radiations than the signal) are discarded. 

Above all, a high significance as 7 to 9a can be achieved to observe WWj productions 
at the 14 TeV LHC, depending on the cuts, e.g., 7.4er, 9.1cr and 8.8cr for cases (A), (B) and 

1 Note the K factor depends on Pr-y according to Ref. [36], where it can be 1.5 for Pt 7 ~ 20 GeV. Here 
we take the K factor for Pr-y ~ 80 GeV. 
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(A) Jet vetoing cut P* (B) P T lower cut (C) Jet reconstructing cut P' 




100 200 3C0 .::>(: :lo:> RIO 900 1000 200 400 600 800 1000 

P T /GeV P T /GeV 



Figure 3. Comparing the differential distributions in leading lepton Pt and Photon Pt, for WWjj 
in the SM, and with QGCs k^ c = ±3 x 10" 5 GcV" 2 

(C), respectively. Note a large portion of WW'y events can come from the QED ISR/FSR 
WW which is not related to QGCs, as shown by the pure_Vs-curves in Fig. 2, however, 
sticking to large Px-y lower cut (~ 80 GeV), one can still get a high significance about 4.5 
sigma and almost totally from pure_Vs contributions. 

5.2 Anomalous WW^fj Couplings 

The W + W _ 7 signal process can be sensitive to anomalous QGCs WW77 and WWZ-y. In 
this paper we are only considering WWjj QGCs for simplicity. As shown in Fig. 3, which 
compares WW^/j in the SM and with k^ c = ±3 x 10 _5 GeV" 2 , the QGCs lead to excesses 
on the hard tails in various kinematic region. One thus can refine the cuts in Sec. 4 to 
enhance the sensitivity to QGCs as following: 

• (a) P T7 > 250 GeV. 

• (b) The leading lepton P T > 200 GeV. 

• (c) \m u - M z \ > 5 GeV. 

After all these selection cuts, the significances are calculated and displayed in Fig. 4 
as functions of the QGCs k^, at the 14 TeV LHC, with an integrated luminosity of 30, 
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Figure 4. Significance with different values of WWj'y anomalous couplings at y/s = 14 TcV LHC 



100 and 200 fb~ . The horizontal dash lines here correspond to the 95% confidence level 
limit. Note here the signal is defined as (QGCs — WW77). 

We also compare our results with the existing constrains given by the OPAL Col- 
laboration [14] and the previous MC expectation limits based on VBF photon exchange 
processes [16], as shown in Table 2. At the 14 TeV LHC, we can set more stringent limit 
than the OPAL, down to 1 x 10~ 5 GeV" 2 level with only 30fb _1 LHC data. Although our 
W + W~ r ) channel seems to set more loose limits on QGCs than the VBF photon exchange 
process [16], however, it has simpler event topology and may be less contaminated by the 
QCD and VBF systematics. 



Couplings 


OPAL 


Photon exchange 


W + W 7 production 


[xlO" 5 GeV- 2 ] 


xlO" 5 


30(200) ft." 1 


30 fb- 1 


100 fb" 1 


200 fb" 1 


«o W /A 2 


[-0.020, 0.020] 


0.26 (0.14) 


[-1.40, 1.42] 


[-0.82, 1.04] 


[-0.66, 0.66] 


af/A 2 


[-0.052, 0.037] 


0.94 (0.52) 


[-2.60, 2.26] 


[-1.86, 1.60] 


[-1.62, 1.24] 



Table 2. 95% C.L. limits for existing LEP2 results [14], the VBF photon exchange process [16] 
and the W + W~~ r y triple boson production process. 



6 Discussion 



In the past, due to the insufficiency of the available center-of-mass energy, measurement on 
triple boson production has never been performed at hadron colliders. With the excellent 
behavior and possible upgration in the near future of the powerful LHC, this kind of 
measurement could become possible in the future. Our MC simulation study shows that 
at the 14 TeV LHC with an integrated luminosity of 100 (30) fb" 1 , one can reach a 
significance over 9 (5) a to observe the SM W + W~j production, and can constrain at the 
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95% CL the anomalous IUVT77 coupling parameters, e.g., Oq / A 2 (see Ref. [14] for their 
definitions), below 1 x 10 _5 GeV~ 2 , respectively. The expected limits are far beyond the 
existing LEP results, although less tighter than the ones from MC studies on VBF photon 
exchange processes, suffer less the QCD and VBF systematics due to cleaner and simpler 
events topology. 
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